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Abstract 


The  purpose  of  this  note  is  to  investigate  some  possibilities  of 
obtaining  in  a  simple  manner  approximate  solutions  of  the  bound  two-electron 
atomic  system.  In  particular,  an  effective  two-body  method  is  studied.  It 
consists  essentially  of  the  assxamption  of  a  simple  functional  form  for  the 
wave  function.  The  resulting  Schrodinger  equation  may  be  solved  exactly, 
yielding  upper  bounds  of  all  the  eigenvalues  of  the  system.  Only  the  lowest 
approximate  eigenvalue  is  near  to  the  exact  result.  Some  methods  of  improve- 
ment are  suggested  and  discussed. 
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1.    Introduction 

In  the  present  note  ve  wish  to  present  some  observations  on  the  non- 
relativistic  two-electron  atomic  system.  We  hope  these  may  prove  xiseful  and 
stimxilating  to  additional  efforts  to  obtain  approximate  eigenfunctions  and 
eigenvalues.  In  particular,  we  examine  the  reduction  of  the  problem  to  an 
"  effective  two-body'  system,  an  idea  that  has  been  used  in  the  corresponding 
three-body  nuclear  system,  the  triton  L  J '  L  J  ^  A  complete  review  of  the  extensive 
treatment  of  the  two- electron  atomic  system,  in  parti cvilar  the  helium  atom,  may 
be  found  in  the  recent  book  of  Bethe  and  Salpeter  l-^-' .  The  variational 
formulation  of  the  problem,  which  is  the  basis  of  the  present  investigation,  is 
reviewed  here. 

Instead  of  the  electronic  radial  coordinates  r  ,  r  and  the  electronic 
interdistance  r^  ,  it  has  been  foimd  convenient  in  the  past  to  use  the  variables 

(1)  s  =  r^  +  r^  }    t  =  r^  -  r^   j     ^  =  ^12 

The  variational  principle  for  the  Schrodinger  eqviation  for  a  two-electron  atom 
in  a  state  of  zero  total  anguleir  momentian  may  be  written  as 

(  ? 

(2)  5  J    {^E^  -  E^ndT    =  0 

where 

(3)  f  m.r  =  r^  r  \u  I  \t  |(s^-t^)u[(||)^  .  (^)=  .  i^)']  .  2s(u^-t2)|i  ^ 

^o  o        ^-u 


*  2t(s2.u2)^  ^  +  (s2.t2-l.Zsu) ^2 


(M 


00  s  u 

l/^dT    =1       ds    I     du  I     dt(s2-t2)ui2 

-u 


2  - 


and 

(5)  ^  =  ^(s,t,u)  . 

The  parameter  Z  is  the  charge  nimibdr  of  the  nucleus,  for  example,  Z  =  2 

for  the  heliimi  atom.   The  units  are  those  of  Hartree:  m  =  1,  e  =  1, 

— ^  =  a  =1.  In  these  units,  the  ground  state  energy  of  the  hydrogen 

me        ^ 

atom  is  (-  p).  The  unit  of  energy  in  electron  volts  is  27.1. 

The  most  accurate  calculations  on  the  ground  state  of  the  heliian  atom 

have  utilized  the  Ritz  variational  form  of  the  energy  with  the  insertion 

of  apprcpriatfe  tirlal  wave  functions  \|f  .  If  the  electronic  interaction 

2  2  2 
(the  term  (e  -t  )\t'  in  equation  (5))  were  absent,  the  exact  ground  state 

wave  function  would  be  the  product  of  two  hydrogenic  ground  state  wave 

functions 

-Zr    -Zr 

(6)  i|r^  =  e   ■^.  e   2  =  e  ^   . 

Consequently,  in  one  of  his  very  early  calculations,  Hylleraas  *■  ■•  chose 
^Ir  as  a  trieJ.  function  with  Z  a  variational  parameter,  and  obteiined  the 


best  value  Z  =  27/l6,  for  which 


(7)       E^  =  -2.85_ 


This  is  already  within  2  %  of  the  ejgperimental  value  of  the  paraheliimi 
( 'S)  ground  state, 

(8)  E     =  -2.90»*  . 
^  '       o  exp 

Many  efforts  by  Hylleraas  and  others  have  been  made  to  improve  this  result 
by  assuming  the  trial  functions  to  be  of  the  form 

(9)  t,  .  e---  E  ^,,.,y  -V 

With  k  and  c„  as  variational  parameters.  The  most  ambitious  calciilation 

jp,m,n 


[51 
along  these  lines  was  performed  recently  by  Kinoshita '■^-'  who  utilized  a  series 

of  58  terms.  His  result  for  the  lowest  energy  value  is 

(io)    ^o  "  '  2.9037225  . 

Such  calculations  require  an  extensive  amount  of  labor.  In  looking  into 
this  problem,  we  were  motivated  by  a  desire  to  obtean  approximate  solutions  of 
the  two-electron  system  with  little  computational  effort. 

2.    The  effective  two-body  method 

In  particular,  we  have  applied  the  technique  of  converting  the  three-body 
problem  to  an  effective  two-body  problem  which  has  been  utilized  in  the 
investigation  of  the  binding  energy  of  the  triton,  the  coirespondlng  three- 
body  nuclear  problem'-  •'''■-'.  The  restriction  of  ^If.  to  be  a  function  of 
only  one  variable  effects  such  a  conversion.  Thus,  if  we  let  ^f.  =  0(s)  only, 
the  variational  principle  (l)  reduces  to 

The  Euler-Lagrange  equation  for  0(s)  is  then 

s 
So  far  «n  this  was  observed  by  Hyileraas,  who  also  noted  that  the  I0w6st 
eigenfunction  and  eigenvalues  are  given  by  equations  (6)  and  (7)^  respectively. 
It  does  not  appear  to  have  been  noticed  previously  that  (12)  may  be  solved  for 
all  its  eigsnfunctions  and  eigenvalues.  Setting  0(s)  =  s"-^'T'(s),  the 
equation  for  F  becomes 

^2   8  ds   ^     s     ^2  ' 
with  ^  =  |^(l6Z-5)  and  X  =  3/2. 


-  k 


This  is  readily  recognized  as  the  hydrogenic  equation  (see  Appendix). 

Consequently,  F  is  to  be  interpreted  as  an  "  effective"  two-body  vave  fimction 

representing  a  hydrogenic  system  for  vhich  the  charge  on  the  nucleus  is  ^ 

and  the  effective  orbital  qtiantum  number  is  X.  The  non-  singular  eigenfunctiona 

of  this  equation  are  sinple  generalizations  of  the  usual  hydrogenic  solutions 

and  may  be  written  as 
_  P 
(Ik)  F^  =  p^  2  L^(^)  ,   p  =  2ks,   k^  =  |E|  , 

with  L  (p)  the  generalized  Lagrange  function  (See  Appendix).  In  particular, 

F  =  cP'^e~  ^'    '     The  associated  eigenvalues,  Just  as  in  the  hydrogen  piroblem, 
o 

are 

(15)       K     = ^ 5  y  n  =  0,1,2 

°     (n+X+1)^ 

It  is  known  that  for  each  n,  the  n-th  largest  eigenvalue  of  a  matrix  which  is 
reduced  by  one  row  and  column  is  an  upper  bound  on  the  n-th  largest 
eigenvalue  of  the  original  matrix.  This  theorem  was  first  applied  in  the 
present  variational  context  by  Hylleraas  and  Undheim"-  ^ ,   and  it  follows 
that  the  above  eigenvalues  are  respective  ttpper  bounds  of  the  true  'S 
eigenvalues  of  the  two-electron  system.  Furthermore,  the  eigenfunctions 
foim  a  ccanplete  orthogonal  set  in  the  variable  s  over  the  domain  (0,oo)» 
This  result  for  the  ground  state  eigenvalue  E  compares  well  with  the 
experimental  values,  as  may  be  seen  in  Table  I.  The  values  of  the  lowest 
ftigenvalueis  for  all  two-electron  atoms  from  Z  =  1  to  Z  =  9  are  given  there. 
Naturally  these  same  resxilts  could  have  been  obtained  by  choosing  as 
a  trial  function  a  simple  exponential.  We  may,  however,  make  the  stronger 
assertion  that  the  exponential  form  is  the  'best'  function  wliLch  is 
exclusively  a  function  of  s.  This  observation  is  independent  of  the  value 
of  Z;  in  fact  we  see  that  explicit  minimization  with  respect  to  Z  has  been 
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avoided.  It  may  also  be  observed  from  Table  I  that  the  discrepancy  from  the 

expertmental  resiilt  is  more  or  less  tonstant  with  Z  so  that,  percentagewise, 

this  error  becomes  more  negligible  as  Z  increases.  Conversely,  the  value 

for  H'  (Z=1)  is  not  even  sufficient  to  bind  the  second  electron  because  E  is 

o 

greater  than  ( -  -p) .  The  veilues  obtained  for  the  second  state  are  uniformly 
poor:  for  all  Z  the  second  state  is  given  as  E  =  (rQ)E  and  this  value  is 
never  below  the  enei^y  level  of  a  single  electro  n  attached  to  a  nucleus  with 
charge  Z.  Thus,  th«»  effective  two-body  approximation  introduces  a  special 
simplicity  -vdiich  may  be  considered  a  good  approximation  to  the  true  wave 
function  only  for  the  groxmd  state.  Nevertheless,  the  'S  solutions  may  be 
useful  not  only  in  providing  a  complete  orthogonal  set  of  functions  in  the 
variable  s,  but  also  in  indicating  the  "best*  form  of  the  eigenJFunctions 
in  their  dependence  on  s. 

5.    Introduction  of  a  free  parameter 

The  ease  with  which  it  is  possible  to  obtain  these  approximate  solutions 

makes  it  appear  desirable  to  search  for  ways  of  improving  the  acctiracy  while 

still  retaining  the  ease  of  solution  of  the  restilting  two-body  Schrbdinger 

equation.  One  way  of  accomplishing  this  is  to  introduce  the  variable  S  =  s-Ktu 

instead  of  s,  where  a  is  an  arbitrary  parameter.  The  volume  element  now 

becomes 

OD     S/l+a    u 

(16)     I  dr  =  I   dS  I      du  I  dt[(s-au)^-t^Ju  with  a  >  -1  . 

-^       o     o        -u 

Assvmie  now  that  i.  =  0(s)  only.  We  may  again  integrate  over  u  and  t  to  obtain 


(17)    5 


00 
f 


dS  , 
o 


a(,)  I,  s5(M)^  .  v(a,Z)  i%2  sV.  n(a4  E  S^^ 


50  "  MS'    -v-'"'   12    >"   -^""50 
where        a(a)  =  |(8  +  15a  +  lix^  +  5a^  +  a  ) 


\   =  0 


6  - 


1      2 
n(a)  =  ^(8+5a+a  )  , 

Clearly,  setting  a  =  0  reduces  this  equation  to  the  one  previously  considered. 
Purthermore,  the  eigenfunctions  and  eigenvalues  axe  the  same  as  before  provided 
ve  substitute  ^(Z)  ^^W   for  ^  and  E  ^p|  for  E.  The  eigenvalues  are 
therefore  given  as 

(^R\  P  -  [5  (l6Z-^n2  v%,Z) 

^^^^       \  -   -[l5  (2n+5lJ   a(a)n(a)  ' 

The  eigenvalues  etre  again 
(19)       0^  =  e-p/V(p) 

vlth 

p  =  2k  S   and  k^  =  l^l^r  E  I . 
n        n   la^a;  nl 

This  expression  for  E  may  be  minimized  with  respect  to  a.,   because  it  still 

constitutes  a  stationary  expression  for  E  .  For  helium,  the  best  v41ue  of 

n 

a  is 


a  =  -  0.157 

(20)      Jforvhlch    E  =-2.889618 

'  °    -k  (s-0.157u) 


and        0  =  e 

'^o 


o 


Of  course,  this  dirplicates  the  resiJLLt  of  Hylleraas  obtained  vlth  the  trial 
function  t*  =  e      vlth  k  and  c  variational  parameters.  The  general 
remarks  made  previously  (for  the  case  a  =  O)  also  apply  here.  The  optimum  value 
of  a  depends  on  Z.  In  Table  I  may  be  found  the  optimvim  values  of  a  together 
with  the  lowest  eigenvalues  for  all  two  electron  atoms  iip  to  Z  =  9.  As  for 
the  previous  case,  the  error  is  more  or  less  independent  of  Z.  The  least  error 


-  7  - 


is  made  for  Z  =  9:  less  then  l/lO  of  l7o  percentage  error  is  made  for  Z  =  1, 
for  vhich  we  see  that  the  H~  atom  is  barely  boxmd.  However  the  corresponding 
eigenfunction  0  as  given  by  (20)  is  by  no  means  the  best  siniple  trial 

function  for  H".  Oiandrasekhar'-  ■'  obtained  E  =  0.5135  with  a  trial 

-ks 
function  of  the  form  t+  =  e    cosh  ct.  Such  a  function  is  not  readily 

generalized  in  a  simple  way.  The  difficulty  really  stems  from  the  Pauli 

principle:  ^  must  be  even  (for  the  'S  states)  or  odd  (for  the  ''s  states) 

in  the  variable  t.  To  generalize  this  in  a  simple  manner  requires  the  assumption 

(21)      ^^   =  0(s  +  pt)  +  0(s  -  pt)  . 

[91 
Assumption  (21)  is  really  an  exaaiple  of  Fock's  method"-  -■  in  a  slightly 

different  guise  and  experience  shows  that  it  leads  to  an  integro-differential 

equation  for  0.  It  may  be  that  the  two-body  assumption  made  here  is 

sufficiently  good  to  make  it  worth-while  to  solve  this  equation  numericeOly. 

A  more  important  reeison  for  considering  this  equation  is  that  a  similar 

asstmrption  is  necessary  to  obtain  a  good  representation  of  the  excited  states*-  ■' . 

k.         Suggestions  for  improvement 

Can  we  generalize  the  effective  two-body  method  here  introduced  and 

still  retain  the  simplicity  of  solution  of  the  final  equation?  It  appears 

that  the  assumption  i|f.  =  0(s)  G(-  ,  -)  where  G  is  any  prescribed  (square 

integrable)  function  of  the  variables  indicated,  leads  to  essentially  the 

variational  principle  (ll),  with  different  values  of  the  numerical  coefficients 

/s   /u 
for  different  choices  of  G.  This  is  easy  to  see  because    du  1  dt 

takes  on  constant  limits  in  terms  of  the  new  variables.  However,  Kinoshita"-  J 

has  shown  that  an  exact  power  series  solution  of  the  gro\and  state  wave 

function  of  the  heliimi  atom  would  not  contain  terms  of  the  ftjim  (— )  (— ) 

where  n,  m  are  positive  integers.  Therefore,  the  additional  freedom  of 


8  - 


choosing  the  function  G  would  not  appear  to  be  loseful,  even  in  the  variational 
principle.  An  improvement  in  accuracy  natirrally  reqtiires  the  introduction 
into  the  wave  function  of  the  dependence  on  the  variables  u  and  t.  A 
natural  attempt  to  improve  the  results  of  the  two-body  method  is  to  let 

(22)        ^^   =  0(s)  •X(u,t) 

with  0(s)  prescribed,  and  to  solve  the  resulting  Euler  equation  for  X* 
Such  an  assumption  may  be  viewed  as  an  •enforced*  separation  of  variables. 
There  wovjld  result  a  Schrodinger  eqiiation  in  only  the  two  independent 
variables  u  and  t.  One  might  hope  to  get  some  insight  in  this  manner  into 
the  dependence  of  the  time  wave  function  on  these  variables.  It  is  simpler 
to  restrict  %  "to  be  a  function  of  u  or  t  above,  in  -which  case  it  too  becomes 
an  effective  two-body  wave  function.  For  exaa^le,  to  improve  the  ground 
state  wave  function  0  (s),  let 


(25) 


*t  =^ 


_  «-k  s/2 


X(u) 


With  distances  normalized  by  k,  there  results,  sifter  integraticn  over  t  and  s 
oo 


{^) 


f 


J 


du.u2  .(u)(M)2.t(,)  xf^.^t-  i[zv(u)-  2^]X^-  \  .(u)  t 


=  0 


where   a(u)  =  e"  (l+u+  u  /5)  , 
b(u)  =  e-"(u  +  u2)/5  , 
v(u)  =  e'^'d  +  u)  . 
Setting  X(^)  =  *>•'  ^(u),  the  Schrodinger  equation  for  F  is 

E 


(25)    4^-^-^ 
,2   u  du 
du 


1  j^  I    \   ^'^tn   v(u)   1  X 
-  5-  +  c(u)  +  -(z  ^  .  ^ 


F  =  0 


00 

subject  to  the  normalization  condition     du 

=(u)   =  k,  i2i%^ 


F^.   u-  =  1. 


The  function 


35 


(l+u+u75)' 

is  fairly  constant  over  the  entire  range  of  u.  This  equation  is  reminiscent 
of  the  sort  of  inter-atomic  potentials  occurring  in  molecules,  but  it  appears 
that  it  would  have  to  be  solved  numerically.  As  in  the  Ritz  method,  the 
optimum  value  of  k  has  to  be  determined  by  variation,  and  the  value  given  by 
{Zh),   say,  is  not  necessaxily  the  best  one.  It  is  perhaps  worth  observing 
that  the  asymptotic  form  of  the  equation  for  F  for  large  values  of  u  is  again 
hydrogenic. 


(26) 


In  a  similar  manner,   the  assxjnrption 
_  k 
>y^  =  e"  2  \  (t) 


res\jlts  in 


(27) 


00  /• 

dt     a(t)(f  )2  -  i  v(Z,t)  t  -  a(t)(^  -  \)  t\  =  0 

o        L 


where  a(t)   =  (5  +  t  )   cosh  t  -  5*  sinh  t     , 

v(Z,t)    =  (|  -  Z)cosh  t  +  (Z  -  l)t  sinh  t   . 

In  deriving  this,  we  have  made  use  of  the  fact  that  X  ("t)  Tsast  have  a 
prescribed  synmetry  (even-ness  for   'S  states)   in  the  variable  t.     letting 
X  (t)   =  ta'-'-'^(t),   there  results 


,^a\  d^  ^  2  dP  ^ 

(28)  —  -^  t  dt  ^ 


E         1  ^     /.  X        1  v(Z,t)" 


F  =  0 


dt^       "  Ik  ^    '  J  00 

with  c(t)  =  \{—)^-   I  —,   a»<i  subject  to  the  boundedness  of  I   dftT^. 

o 


-  10 


It  wovild  be  interesting  to  solve  numerically  either  (25)  or  (28)  in 
order  to  compare  vith  some  of  the  known  variational  calculations.  No  doubt, 
it  is  too  much  to  expect  to  improve  upon  Kinoshita's  woiic  in  thig  manner 
because  even  thoijgh  we  have  (in  comparison  with  assumption  (9))  retained 
all  powers  of  u  and  t  (or  both,  if  assumption  (22)  is  used),  we  have  neglected 
additional  powers  of  s.  However,  con?>ared  with  two-electron  oidculations 
which  have  not  been  so  extensive,  it  is  quite  possible  that  this  procediure 
may  be  the  simplest  way  to  get  very  good  approximate  results.  If  the  n 
eigenstate  is  desired,  than  of  covirse  0  (s)  in  (22)  should  be  specified, 
and  the  resulting  equation  for  %  should  be  solved  for  its  nth  eigenstate. 
Finally,  we  mention  that  we  may  combine  the  method  of  Section  2  with  the 
method  above  (restricting  X  "to  be  a  function  of  t  above)  to  obtain 
dependence  of  the  trial  wave  function  on  all  three  variables. 


-  11  - 


Appendix 


51ie  radial  equation  for  a  hydrogenic  atom  with  charge  nmnber  Z  and 
orbital  quantum  number  X  is 


d^  ^  2  dE  ^ 
dr^   ^^ 


2E 


+  2Z  _  X(X+1) 


R  =  0  . 


The  eigenfunctions  for  the  bound  states  are 

P 


where 


R„  =  p\  ^L^(P) 


p  =  2kr 


2E  =  -  k 
and  L  is  the  generalized  Laguerre  function*-  -' 

0=0 
They  satisfy  the  relationship 

_,  _X/    \  p  -X  d       /    -p  n+Xv 

nl  Ljp)   =  e^p       — J  (e  ^p       ). 

dp 
The  first  few  generalized  Laguerre  functions  for  fixed  X  are 


L^(P)    =  1 

L^(p)    =  X  +  1  -    p 


LgCp)    =  |[(X+l)(X+2)    -   2(X+2)p  +  p2]      . 


The  associated  eigenvalues  are  E 


1  Z^ 

2 


.2  • 


(n+X+1)' 
In  the  notation  used  here,  n  is  the  radial  quantum  number  and  not,  as  is 

frequent  T;i6age,  the  total  quantim  number  (n+X+l). 
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